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ABSTRACT 

Explicit  algebraic  stress  models  that  are  valid  for  three-dimensional  turbulent  flows  in 
non-inertial  frames  are  systematically  derived  from  a  hierarchy  of  second-order  closure  mod¬ 
els.  This  represents  a  generalization  of  the  model  derived  by  Pope  [ J .  Fluid  Mech.  72,  331 
(1975)]  who  based  his  analysis  on  the  Launder,  Reece  and  Rodi  model  restricted  to  two- 
dimensional  turbulent  flows  in  an  inertial  frame.  The  relationship  between  the  new  models 
and  traditional  algebraic  stress  models  -  as  well  as  anistropic  eddy  viscosity  models  -  is 
theoretically  established.  The  need  for  regularization  is  demonstrated  in  an  effort  to  explain 
why  traditional  algebraic  stress  models  have  failed  in  complex  flows.  It  is  also  shown  that 
these  explicit  algebraic  stress  models  can  shed  new  light  on  what  second-order  closure  models 
predict  for  the  equilibrium  states  of  homogeneous  turbulent  flows  and  can  serve  as  a  useful 
alternative  in  practical  computations. 
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tract  Nos.  NAS1-18605  and  NAS1-19480  while  the  second  author  was  in  residence  at  the  Institute  for 
Computer  Applications  in  Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton, 
VA  23681. 


1.  INTRODUCTION 


Turbulent  flows  contain  length  and  time  scales  that  can  change  dramatically  from  one 
flow  configuration  to  the  next.  Consequently,  two-equation  turbulence  models  -  wherein 
transport  equations  are  solved  for  two  turbulent  fields  that  are  directly  related  to  the  length 
and  time  scales  -  represent  the  simplest  level  of  Reynolds  stress  closure  that  can  be  geometry 
independent.  This  explains  the  current  popularity  of  the  K  —  e  model  which,  unlike  the  older 
mixing  length  models,  does  not  require  the  specification  of  empirical  length  or  time  scales 
that  must  be  adjusted  in  an  ad  hoc  fashion  from  one  flow  to  the  next  (see  Launder  and 
Spalding  1974).  However,  despite  this  positive  feature,  the  K  —  £  model  shares  a  common 
deficiency  with  the  older  mixing  length  models:  it  is  based  on  the  Boussinesq  eddy  viscosity 
hypothesis.  It  is  well-known  that  eddy  viscosity  models  are  unable  to  properly  describe 
turbulent  flows  with  body  force  effects  arising  from  curvature  or  a  system  rotation.  Likewise, 
their  inaccurate  prediction  of  normal  Reynolds  stress  differences  make  eddy  viscosity  models 
incapable  of  describing  secondary  flows  in  non-circular  ducts  (c.f.  Lumley  1978,  Launder  1990 
and  Speziale  1991).  This  has  led  to  the  development  of  nonlinear  stress-strain  relationships 
in  turbulence  modeling  that  transcend  the  Boussinesq  hypothesis  -  an  area  of  research  that 
has  received  varying  degrees  of  attention  during  the  past  three  decades. 

Early  work  on  the  development  of  turbulence  models  with  a  nonlinear  stress-strain  rela¬ 
tion  tended  to  be  empirical  in  nature,  relying  on  analogies  with  non-Newtonian  flows  (see 
Rivlin  1957,  Lumley  1970  and  Saffman  1977).  The  similarities  between  the  laminar  flow  of  a 
non-Newtonian  fluid  and  the  mean  turbulent  flow  of  a  Newtonian  fluid  have  long  been  rec¬ 
ognized  (see  Hinze  1975  for  an  interesting  discussion  of  this  point).  In  these  older  empirical 
models,  which  typically  were  obtained  by  simple  tensor  invariance  arguments,  the  Reynolds 
stresses  were  taken  to  be  nonlinear  polynomial  functions  of  the  mean  velocity  gradients.  In 
recent  years,  nonlinear  Reynolds  stress  models  of  this  type  have  been  obtained  within  the 
context  of  two-equation  turbulence  modeling  by  more  formal  expansion  techniques  incorpo¬ 
rating,  for  example,  the  Direct  Interaction  Approximation  (DIA)  and  the  Renormalization 
Group  (RNG)  (see  Yoshizawa  1984,  Speziale  1987,  Rubinstein  and  Barton  1990,  and  Yakhot 
et  al.  1992).  These  models,  which  are  characterized  by  an  explicit  relationship  between  the 
Reynolds  stress  tensor  and  the  mean  velocity  gradients  (and  possibly  their  time  derivatives) 
have  come  to  be  referred  to  as  “anisotropic  eddy  viscosity  models.” 

During  the  1970’s,  second-order  closure  models  became  popular  in  which  closure  was 
achieved  based  on  the  Reynolds  stress  transport  equation  where  history  and  nonlocal  effects 
are  accounted  for.  By  means  of  an  equilibrium  hypothesis,  in  which  the  Reynolds  stress 
convection  and  transport  terms  were  neglected,  algebraic  stress  models  were  obtained  from 
second-order  closures  (see  Rodi  1976).  In  these  models,  the  Reynolds  stresses  were  related 
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implicitly  to  the  mean  velocity  gradients.  This  provided  the  first  formal  means,  based  on 
a  higher-order  closure,  to  justify  the  extension  of  the  Boussinesq  hypothesis  to  incorporate 
nonlinearities  in  the  mean  velocity  gradients.  However,  this  type  of  algebraic  stress  model 
is  cumbersome  to  implement  in  complex  flows  since  the  stress-strain  relation  is  not  explicit; 
numerical  stiffness  problems  can  result  from  the  need  for  successive  matrix  inversions  at  each 
iteration. 

In  an  interesting  paper  that  has  to  a  large  extent  gone  unnoticed,  Pope  (1975)  developed 
a  nonlinear  Reynolds  stress  model  by  invoking  the  same  equilibrium  hypothesis  as  Rodi 
(1976).  However,  Pope  (1975)  actually  presented  a  methodology  for  obtaining  an  explicit 
relation  for  the  Reynolds  stress  tensor  from  the  implicit  algebraic  equations  that  Rodi  (1976) 
obtained  from  the  Launder,  Reece  and  Rodi  (1975)  model.  This  methodology,  which  leads 
to  explicit  algebraic  stress  models,  is  based  on  the  use  of  integrity  bases  from  linear  algebra. 
Due  to  the  complexity  of  the  algebra,  Pope  (1975)  was  only  able  to  obtain  a  solution  for 
two-dimensional  turbulent  flows  which  he  accomplished  using  the  Launder,  Reece  and  Rodi 
model. 

The  main  purpose  of  the  present  paper  is  to  extend  the  results  of  Pope  (1975)  to  three- 
dimensional  turbulent  flows  in  non-inertial  frames  starting  from  a  more  general  hierarchy  of 
second-order  closure  models  -  a  task  that  is  now  feasible  computationally  via  symbolic  manip¬ 
ulation.  The  relationship  between  these  new  explicit  algebraic  stress  models  and  anisotropic 
eddy  viscosity  models,  as  well  as  the  older  nonlinear  Reynolds  stress  models,  will  be  estab¬ 
lished  in  a  systematic  fashion.  There  appears  to  be  considerable  confusion  in  the  turbulence 
literature  concerning  the  relationship  between  these  various  types  of  nonlinear  Reynolds 
stress  models  that  needs  to  be  clarified.  It  will  be  shown  that  the  new  explicit  algebraic 
stress  models  represent  the  equilibrium  Reynolds  stress  anisotropies  predicted  by  second- 
order  closures  in  homogeneous  turbulent  flows.  Consequently,  beyond  their  potential  use  in 
practical  turbulence  calculations,  the  explicit  algebraic  stress  models  to  be  derived  herein 
can  be  used  to  explore  the  predictive  capabilities  of  a  hierarchy  of  second-order  closures  in 
homogeneous  turbulence.  These  issues  will  be  discussed  in  more  detail  in  the  sections  to 
follow. 

2.  THEORETICAL  BACKGROUND 

The  incompressible  turbulent  flow  of  a  viscous  fluid  will  be  considered  where  the  velocity 
v  and  kinematic  pressure  P  are  decomposed  into  ensemble  mean  and  fluctuating  parts  as 
follows: 

v  =  v  T  u,  P  =  P  +  p.  (1) 
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The  Reynolds  stress  tensor  r,j  =  U{Uj  is  a  solution  of  the  transport  equation  (c.f.  Hinze 
1975): 

Dt' '  •  $V' 

—  Ufc  Qx  Tjfc  Qx  2f ^‘mi^mkjTik  "I"  ^mki^jk )  "h  ^ij  "1"  A?  (2) 

which  is  valid  in  an  arbitrary  non-inertial  reference  frame  that  can  undergo  a  rotation  with 
angular  velocity  fl,  relative  to  an  inertial  framing.  In  (2),  D/ Dt  =  d/dt  -f  v  •  V  is  the 
convective  derivative,  eijk  is  the  permutation  tensor,  v  is  the  kinematic  viscosity,  and 


<b 
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(5) 


are,  respectively,  the  pressure-strain  correlation,  the  dissipation  rate  tensor,  and  the  turbu¬ 
lent  transport  term.  Turbulence  models  based  on  the  Reynolds  stress  transport  equation  (2) 
are  referred  to  as  “second-order  closures”  since  (2)  is  obtained  by  taking  a  second  moment 
of  the  fluctuating  Navier-Stokes  equation. 

Homogeneous  turbulent  flows  in  equilibrium,  as  well  as  regions  of  inhomogeneous  turbu¬ 
lent  flows  where  there  is  a  production-equals-dissipation  equilibrium,  satisfy  the  constraints: 


where 
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(7) 

(8) 


is  the  anisotropy  tensor  (K  =  fuiui  is  the  turbulent  kinetic  energy).  In  physical  terms, 
this  is  an  equilibrium  for  which  convective  and  transport  effects  can  be  neglected;  it  is  the 
basic  equilibrium  hypothesis  used  in  the  derivation  of  algebraic  stress  models.  Although 
it  constitutes  an  idealization,  it  is  comforting  to  know  that  this  equilibrium  hypothesis 
is  achievable  in  interesting  physical  limits  that  include  homogeneous  shear  flow  and  the 
logarithmic  region  of  an  equilibrium  turbulent  boundary  layer. 

It  follows  from  (6)  that 
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□  □ 


(10) 


and,  hence,  by  making  use  of  the  contraction  of  (2)  and  (7),  we  obtain 

=  (v  -  £)hi 

Dt  [  ’  K 

where  V  =  —Tijdvi/dxj  is  the  turbulence  production  and  £  =  |£„  is  the  scalar  turbulent 
dissipation  rate.  The  substitution  of  (7)  and  (10)  into  (2)  yields  the  following  equilibrium 
form  of  the  Reynolds  stress  transport  equation: 

t  ■  c)\j  dv 

^ )  wr  Tik  TZ  Tjfc  1"  2,£lm(cmkjTik  T  GmkiTjk)-  (11) 

A  OXk  OXk 

The  dissipation  rate  tensor  can  be  split  into  isotropic  and  deviatoric  parts  as  follows 

2 

£\j  =  d"  D^ij-  (12) 

By  making  use  of  (8)  and  (12),  we  can  rearrange  (11)  into  the  alternative  form 


(V-e)btJ  =  --KSij  —  K(bikSjk  +  bjkSik 

^mn^mn^ij)  ^  \bik{^jk  "T  2em^jfim) 


-\-bjk(u?ik  +  2emk,f}n)]  +  -  n,j 

where 

•g  _  l  fdvj  dvj  \  _  1  /  dvj  dvA 

'3  2  \dxj  dxi)  ’  U'3  2  dxi) 

II,j  =  D&ij-  (15) 

In  all  of  the  commonly  used  second-order  closure  models,  Tly  is  modeled  in  the  general  form 
(see  Lumley  1978,  Reynolds  1987  and  Speziale  1991): 

n  a  =  eAijih)  +  KMijki{h)j±.  (16) 

For  non-inertial  frames,  dvk/dxi  is  replaced  with  dvk/dxi  +  in  (16).  The  substitu¬ 

tion  of  (16)  (with  its  non-inertial  correction)  into  (13)  yields  a  closed  system  of  algebraic 
equations  for  the  determination  of  the  Reynolds  stress  anisotropy  in  terms  of  the  mean  ve¬ 
locity  gradients.  This  constitutes  the  general  form  of  algebraic  stress  models;  the  algebraic 
stress  model  of  Rodi  (1976)  is  obtained  when  the  Launder,  Reece  and  Rodi  model  for  IltJ 
is  introduced  into  (13)  and  non-inertial  effects  are  neglected.  It  is  clear  that  these  algebraic 
stress  models  are  implicit  in  nature  since  the  Reynolds  stress  tensor  appears  on  both  sides 
of  the  equation.  In  the  next  section  we  will  examine  how  explicit  relations  can  be  obtained 
from  (13). 
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3.  EXPLICIT  ALGEBRAIC  STRESS  MODELS 


When  models  for  II, j  that  are  tensorially  linear  in  the  anisotropy  tensor  6,j  are  considered, 
it  is  possible  to  obtain  an  explicit  expression  for  b{j  in  terms  of  the  mean  velocity  gradients. 
This  can  be  accomplished  by  techniques  from  linear  algebra  that  make  use  of  integrity  bases 
(namely,  polynomial  representations  for  isotropic  tensor  functions).  Pope  (1975)  was  the  first 
to  show  this  by  basing  his  analysis  on  the  Launder,  Reece  and  Rodi  model  which  is  the  most 
commonly  used  model  for  II, j  that  is  a  linear  function  of  6,j.  We  will  extend  the  analysis  of 
Pope  (1975)  to  a  more  general  hierarchy  of  linear  models  for  II, j  where  non-inertial  effects  are 
accounted  for.  Furthermore,  our  analysis  will  be  conducted  for  three-dimensional  turbulent 
flows.  Pope  (1975)  restricted  his  analysis  to  two-dimensional  turbulent  flows  because  of  the 
severe  complexity  of  the  algebra  -  an  obstacle  that  we  will  overcome  by  using  Mathematica™ 
(Wolfram  1988). 

It  can  be  shown  that  the  most  general  form  of  (16)  that  is  tensorially  Imear  in  the 
anisotropy  tensor  6,j  is  given  by  (c.f.  Reynolds  1987  and  Speziale  1991): 


ILj 


—Cxsbjj  +  C2I\  Sij  +  C3h  (^bikSjk  +  bjkSik  —  —b  mn  Smnbij^j 

+C4K(bikWjk  +  bjkWik) 


(17) 


where 

VF ,'j  —  j  T  £mji (18) 

is  the  absolute  vorticity  tensor  (namely,  the  vorticity  tensor  u;,j  relative  to  an  inertial  frame). 
The  coefficients  C\  —  C4  can  be  functions  of  the  invariants  of  bi:  and  can  depend  on  V/e  (of 
course,  (16)  requires  that,  at  most,  C\  is  a  linear  function  of  V/e).  The  model  of  Launder, 
Reece  and  Rodi  (1975),  the  model  of  Gibson  and  Launder  (1978)  and  the  linearized  model 
of  Speziale,  Sarkar  and  Gatski  (1991)  -  where  the  quadratic  part  of  the  slow  pressure-strain 
correlation  is  neglected  -  are  special  cases  of  (17).  For  these  models,  we  have  the  following 
coefficients: 


Launder,  Reece  and  Rodi  Model 


Cx  =  3.0,  C2  =  0.8,  C3  =  1.75,  C4  = 


Gibson  and  Launder  Model 

Cx  =  3.6,  C2  =  0.8,  C3  =  1.2,  C4  =  1.2  (20) 
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Speziale,  Sarkar  and  Gatski  Model 


Cx  =  3.4  +  1.8P/e,  C2  =  0.8  -  1.3 III'2,  C3  =  1.25,  C4  =  0.40 


where  I h  =  bijbij.  It  should  be  noted  that  for  sufficiently  small  anisotropies,  more  compli¬ 
cated  models  can  be  approximated  by  the  form  (17). 

The  direct  substitution  of  (17)  into  (13)  yields  the  equation: 

bij  —  2^t  {  {p2  ~  3)  Sij  +  (C3  —  2)  (bikSjk  +  bjkSik  —  - bmnSmn6ij ^ 

C  -  4  C  4  (22) 

~b(^4  2)  (^'jk  T  ^  ^fnkj  4*  bjk  -J-  —  j  ^ 


where 


If  we  introduce  the  dimensionless,  rescaled  variables: 

S'i,  =  \gr(2  -  C3 )Sti 


K  =  \sr(2  -  C,)  [u>0  +  (^-|)  emi,nm 


then  (19)  reduces  to  the  simpler  form 


»«  =  -S'j  -  (%Sj*  +  i-ts-k  -  ?6:,s;,<„)  +  b-kw;,  +  b-hw;t. 


In  matrix  form,  (28)  can  be  written  as 


b*  =  -S*  -  (b*S*  +  S*b*  -  ^{b*S*}I)  +  b*W*  -  W*b*  (29) 

where  {•}  denotes  the  trace  and  I  denotes  the  unit  tensor.  This  is  a  linear  algebraic  equation 
for  the  determination  of  b*  in  terms  of  S*  and  W*;  the  solution  to  (29)  is  of  the  general 
form 

b*  =  f(S*,W*).  (30) 

Form  invariance  under  an  orthogonal  coordinate  transformation  (with  rotation  tensor  Q) 
requires  that 

Qf  (S*,  W*)Qt  =  f(QS*Qr,  QW*Qt)  (31 ) 
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i.e.,  that  f  be  an  isotropic  tensor  function  of  its  arguments.  It  can  be  shown  that  the 
satisfaction  of  (31)  requires  that 

b*  =  <?(A) T(A)  (32) 

A 

where  is  the  integrity  basis  for  functions  of  a  symmetric  and  antisymmetric  tensor  and 
G<A>  are  scalar  functions  of  the  irreducible  invariants  of  S*  and  W*.  For  the  case  under 
consideration,  the  integrity  basis  is  given  by  (see  Spencer  1971  and  Pope  1975): 


=  s*w*  -  w*s* 

=  S*2  -  i{S*2}I 

=  w*2  -  i{w*2}i 
=  w*s*2  -  S*2W* 


T(6)  =  W.2S.  +  S.W.2  _  |{S*W*2}I 

t(7)  =  w*s*w*2  -  W*2S*W* 

T{8)  =  s*W*S*2  -  S*2W*S* 

T(9)  =  W*2S*2  +  S*2W*2  -  §{S*2W*2}I 
T<10)  =  W*S*2W*2  -  W*2S*2W* 


The  irreducible  invariants  of  S*  and  W*  are: 

m  =  {S*2}  (34) 

V2  =  {W*2}  (35) 

r?3  =  {S*3}  (36) 

t,4  =  {S*W*2}  (37) 

7,5  =  {S*2W*2}.  (38) 

Hence,  we  have 

G(A)  =  G(a)(t7i,  772,  773, 774, 775),  A  =  1 , 2, ...  10.  (39) 

Pope  (1975)  showed  that  for  two-dimensional  flows,  only  and  as  well  as  7,7 

and  7,2,  are  independent.  Consequently,  the  calculations  become  much  simpler.  Here,  we 
will  obtain  the  three-dimensional  solution  using  the  same  methodology  that  Pope  used 
The  direct  substitution  of  (32)  into  (28)  yields 

£  G(A)T(A)  =  -  £  <*nT(A)  -  £  G(a)[T(a)S*  +  S*T(A)  -  -{T(a)S*}I 


_<p(A)w*  -j.  W*T(A>] 


where  we  have  made  use  of  the  fact  that  =  S*.  Since  T^A^  is  an  integrity  basis  we  know 
that 

T(A)S*  +  S*T(a)  -  -{T(A)S*}I  =  tfA-yT^  (41) 

3  7 
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(42) 


T(A)W*  -  W*T(A)  =  £  J,v>T('y) 

•7 

(i.e.,  any  polynomial  in  S*  and  W*  can  be  expanded  in  the  integrity  bases;  see  Spencer  1971 
and  Pope  1975).  The  10  x  10  matrices  B\y  and  J\y  are  obtained  by  successive  applications 
of  the  Cayley- Hamilton  Theorem  (the  elements  of  these  matrices  are  given  in  Appendix  A). 
By  substituting  (41)  and  (42)  into  (40),  we  obtain  an  equation  from  which  6’*A*  can  be 
determined: 

g{X)  =  -six  -  £  c^H,x  +  £  a^J,x.  (43) 

7  7 

This  is  a  10  x  10  linear  system  of  equations  for  the  determination  of  the  (7*  9  which  can  be 
written  in  the  matrix  form 

AG  =  B  (44) 

where  the  components  of  A,  G  and  B  are  given  by 

A^\  —  +  J\~,  (45) 

G\  =  6,(A)  (46) 

Bx  =  6ia  (47) 

(see  Appendix  A).  The  solution  to  (44)  is  given  by 

a('>  =  a;;.  (is) 

In  order  to  obtain  a  closed  form  expression  for  G’O,  it  is  necessary  to  analytically  invert  the 
matrix  A.  The  manual  inversion  of  matrix  A  is  not  feasible  in  light  of  its  highly  complex 
structure;  however,  such  a  mathematical  computation  can  now  be  done  symbolically  by 
Mathematica™  (Wolfram  1988).  The  resulting  inverse  matrix  A-1  is  quite  complex,  but  as 
noted  in  (48),  only  the  first  column  of  its  elements  are  needed.  The  resulting  solution  for 
<7<A>  is: 


G^  =  (6  -  37,,  -  21r,2  -  2t,3  +  307,4 ) /D 

(7^21  =  —  (3  -f  3t,i  —  67,2  +  27,3  +  67,4)  /  D 
6,(3)  =  (6  -  37,,  -  127,2  -  2 7,3  —  67,4 )  /D 
6,(4)  =  -3  (37,,  -f  2t,3  +  67,4)  /D 
<7<5>  -  9/D 


r;<6>  =  9/d 
G<7)  =  9/D 
6,(8)  =  9/D 
G(9}  =  18/D 
(7<10)  =  0 


where  the  denominator  D  is  given  by 


7  15  2 

D  =  3-  -7,,  +T)1  -  —  7,2  -  87, ,7,2  +37,2  -  7,3+  “7,,  7, 3 


-27,27,3  +  21  7,4  4-  247,5  +  2?,,  7,4  -  6?,2?,4 
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Unfortunately,  this  three-dimensional  form  is  significantly  more  complex  than  the  two- 
dimensional  form  derived  by  Pope  (1975).  Nevertheless,  the  result  presented  in  (49)  can  be 
shown  to  reduce  to  Pope’s  form  in  the  two-dimensional  limit  where  the  invariants  7/3  and  7/4 
are  zero  and  r,  becomes  equal  to  \ii\g-2-  Using  this  ieduction  and  the  redundancies  in  the 
integrity  L.^is  (T*4*  ■  •  -T*10^  are  linear  combinations  of  T^,T^  and  T®3*;  see  Pope  1975). 
the  resulting  expression  for  in  the  two-dimensional  limit  becomes 


6,(l)  =  - 

r;(2)  =  - 

r;(3) 


3 


3  -  2v\  - 

:] 

3  -  2m  -  67/2 

6 


(.51) 


3  -  2 r/i  -  6t/2 


The  anisotropy  tensor  b*  is  then  given  by 

3 


b*  =  - 


3  -  27/1  -  67/2 


S’  -I-  (S’W*  -  W’S’)  -  2  (smi  -  ^{S’2}l)j 


m 


which  can  be  shown  to  be  identical  to  the  expression  in  Pope  (1975)  after  using  (25)-(27) 
along  with  the  constants  (19)  of  the  Launder,  Reece  and  Rodi  model.  It  should  be  noted 
that  Pope  used  an  alternative  form  for  the  third  basis  in  (52)  -  constructed  from  the  two 
and  three  dimensional  Kronecker  delta  -  that  is  directly  proportional  to  our  T*3). 

It  is  now  clear  that  the  model  derived  by  Pope  (1975)  (which  is  a  special  case  of  (52)) 
actually  constitutes  the  explicit  solution  of  the  algebraic  stress  model  of  Rodi  forb  in  terms 
of  S  and  u ;.  Whereas  Pope’s  analysis  is  only  valid  for  two-dimensional  flows  in  an  inertial 
frame,  we  have  now  succeeded  -  by  means  of  the  representation  theorem  embodied  in  (32) 
-  to  extend  his  result  to  three-dimensional  flows  in  non-inertia]  frames. 

Some  comments  are  in  order  concerning  how  the  three-dimensional  solution  obtained 
herein  compares  with  the  result  derived  recently  by  Taulbee  (1992).  The  three-dimensional 
result  of  Taulbee  only  contains  the  bases  T*^T^,T*4*,T^,T*')  where  the  coefficients 
6,(A),  for  the  most  part,  only  depend  on  the  invariant  7/2.  This  simplification  results  from  the 
neglect  of  the  second  strain-dependent  term  in  the  brackets  on  the  r.h.s.  of  (22)  (namely, 
it  is  obtained  by  setting  (  3  =  2).  We  feel  that  this  is  a  questionable  approximation.  Our 
research  on  homogeneous  shear  flow  indicates  that  the  coefficient  (  3  —  2  is  nonzero  and  must 
be  nearly  half  as  large  as  the  coefficient  C4  —  2  in  order  to  obtain  a  good  description  of  shear 
flows.  Furthermore,  by  neglecting  this  term  in  (22),  the  coefficients  of  the  solution  become 
virtually  independent  of  the  strain  rate  StJ  and  the  model  collapses  to  the  standard  K  —  £ 
model  in  the  absence  of  rotational  strains.  This  stands  in  serious  contradiction  of  the  exact 
two-dimensional  result  (52)  as  well  as  our  new  three-dimensional  result. 
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4.  COMPARISON  WITH  OTHER  NONLINEAR  REYNOLDS  STRESS 
MODELS 


The  explicit  algebraic  stress  model  that  we  derived  in  Eqs.  (32)  and  (49)  formally  con¬ 
stitutes  the  equilibrium  anisotropy  tensor  predicted  by  a  hierarchy  of  second-order  closure 
models  for  homogeneous  turbulence.  These  second-order  closures  are  typically  solved  in  con¬ 
junction  with  a  modeled  transport  for  the  turbulent  dissipation  rate  that  is  of  the  general 
form 


£l-C  -V-C  -  +  C  —  (~t  — 

ru  ~  C  fl  r.'''  L'  +  ('e  I  £Tt}Qx. 


Dt  ~e,A"  "‘K  'edxt'  -  l3»-  1  (53) 

where  Ce\,Ce-2  and  Ce  are  either  constants,  or  functions  of  the  invariants  of  b,  that  vary 
from  one  model  to  the  next.  For  homogeneous  turbulence,  the  transport  term  on  the  r.h.s. 
of  (53)  vanishes  and  the  turbulent  kinetic  energy  is  a  solution  of  the  transport  equation 


DK 

Dt 


=  V-£. 


(54) 


From  (53)-(54),  it  follows  that  V/e  achieves  the  following  equilibrium  value  (see  Speziale 
1991) 

-  1 


v_ 

e 


a 


(«) 


that  will  be  used  to  determine  g  in  homogeneous  turbulence.  Hence,  the  explicit  algebraic 
stress  model  derived  herein  can  be  utilized  to  systematically  explore  the  equilibrium  states 
predicted  by  a  hierarchy  of  second-order  closures  in  homogeneous  turbulent  flows. 

As  discussed  earlier,  the  explicit  algebraic  stress  model  given  by  Eqs.  (32)  and  (49) 
represents  an  exact  solution  to  the  general  implicit  algebraic  stress  model  (22)  that  includes 
the  model  of  Rodi  (1976)  as  a  special  case.  Consequently,  by  a  direct  analysis  of  this  explicit 
model,  new  insights  can  be  gained  concerning  the  problems  that  the  traditional  algebraic 
stress  models  have  had  in  the  calculation  of  complex  turbulent  flows.  The  denominator  of 
the  coefficients  G ^  given  in  (49)  contains  a  sum  of  positive  and  negative  terms  which  has 
the  potential  to  become  zero  rendering  singular  behavior.  For  equilibrium  homogeneous 
turbulence  this  will  not  happen  since  rft...rjs  are  constrained  by  the  transport  equations 
for  TtJ  and  £  to  yield  a  well  behaved  solution  (see  Speziale  and  Mac  Giolla  Mhuiris  1989). 
However,  when  an  algebraic  stress  model  is  applied  to  complex  non-equilibrium  turbulent 
flows,  where  the  underlying  assumptions  invoked  in  its  derivation  are  no  longer  valid,  singular 
behavior  can  result.  This  becomes  a  distinct  possibility  in  turbulent  flows  that  have  localized 
strain  rates  that  are  large.  When  the  traditional  algebraic  stress  models  of  the  form  (22)  are 
applied  to  turbulent  flows  with  large  strain  rates  where  D  can  approach  zero,  we  would  expect 
the  model  to  be  ill-behaved.  This  problem  has  been  experienced  in  the  practical  calculation  of 
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complex  turbulent  flows  where  the  Rod!  algebraic  stress  model  has  at  times  failed  to  converge 
when  solved  iteratively  (Demuren  and  Rodi  1984  and  Demuren,  private  communication). 
Hence,  it  is  clear  that  there  is  the  need  to  regularize  these  explicit  algebraic  stress  models. 
Such  an  approach  has  been  used  recently  in  the  kinetic  theory  with  considerable  success  (see 
Rosenau  1989). 

The  model  derived  herein  can  be  regularized  by  a  Pade  approximation  (this  approach 
was  recently  used  by  Yakhot  et  al.  1992  in  the  derivation  of  their  RNG  based  dissipation 
rate  equation).  The  main  idea  behind  this  approach  can  be  easily  illustrated  for  the  two- 
dimensional  model  given  in  (52).  For  this  case,  the  coefficient  on  the  r.h  s.  of  (52)  can  be 


rewritten  as 


where 


3  _  3 

3  —  2??]  —  6772  3  —  2  t?2  +  6£2 


n  =  (S-tf,)'/1,  c  =  (57) 

It  is  clear  from  (56)  that,  for  sufficiently  large  strain  rates  77,  singularities  can  occur;  on  the 
other  hand,  the  rotational  strains  do  not  cause  any  problems.  In  an  equilibrium  homogeneous 
turbulent  flow,  rj  is  typically  less  than  one.  For  example,  in  homogeneous  shear  flow  (with 
shear  rate  .5),  the  parameter  SK/e  achieves  an  equilibrium  value  of  approximately  5.  For 
the  Launder,  Reece  and  Rodi  model,  this  corresponds  to  77  «  0.2  (see  Speziale  and  Mac 
Giolla  Mhuiris  1989).  Hence,  we  want  to  replace  (56)  with  a  regular  function  that,  for  rj 
sufficiently  less  than  one,  is  approximately  the  same.  This  can  be  accomplished  by  a  Pade 
approximation.  For  example,  to  the  first  order  we  take 


7  ‘  1+7?2' 

Substituting  (58)  into  (56)  then  yields  the  expression 

3  *  3JIW) 

3-2771-67/2  3  +  t/2  +  6C2t?2  -P  6C2 

which  is  regular.  Hence,  we  obtain  a  regularized  model  that  is  well-behaved  for  all  strain 
rates  and  constitutes  a  good  approximation  to  the  original  model  within  the  equilibrium 
range  that  it  formally  applies.  For  practical  applications,  the  regularized  model 

b‘  =  -HjwW  [*' +  (s'w'  -  W‘S')  -  2  (S‘2  -  5{S‘2)I)]  (60) 

should  be  used  instead  of  the  form  (52)  first  derived  by  Pope  (1975)  for  the  special  case 
of  the  Launder,  Reece  and  Rodi  model.  For  the  three-dimensional  case,  a  more  systematic 
Pade  approximation  to  (50)  must  be  made.  This  leads  to  rather  cumbersome  expressions 
that  require  further  simplification  and  are  currently  under  investigation. 


s*  +  (S*W*  -  W*S*)  -  2  (V2  -  j{S*2}l)] 
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We  will  now  show  the  relationship  between  the  explicit  algebraic  stress  model  derived 
herein  and  existing  anisotropic  eddy  viscosity  models.  It  is  clear  that  (32),  with  the  coeffi¬ 
cients  (49),  is  of  the  general  mathematical  form 


bh  =  Aw 


(61) 


where  ( dvk/dxi  r  =  su  +  W£[  and  Aijki  is  a  fourth  rank  tensor  that  is  a  function  of  S’  and 
W*.  This  constitutes  the  general  form  of  anisotropic  eddy  viscosity  models;  the  standard 
eddy  viscosity  hypothesis  of  Boussinesq  is  recovered  in  the  limit  as  A^ki  becomes  an  isotropic 
tensor. 

Our  explicit  algebraic  stress  model  for  three-dimensional  turbulent  flows  is  tensorially 
quartic  in  the  mean  velocity  gradients  with  coefficients  that  are  ratios  of  polynomials  in  the 
invariants  of  S  and  W.  We  will  make  use  of  the  fact  that 

T,J  =  ^K6tJ  +  2alKb;j  (62) 

where  c*i  =  ( ' C-i  —  /(C3  —  2).  The  normalized  anisotropy  tensor  b*  can  then  be  expanded 

in  a  Taylor  series.  To  the  first  order  in  the  mean  velocity  gradients,  we  have 


bh  =  -Sb 


(63) 


from  the  algebraic  stress  model  (32)  and  (49).  The  direct  substitution  of  (63)  into  (62),  after 
making  use  of  (25),  yields 

=  (64) 

where  C*  —  \g  (f  —  C2)  •  This  is  the  eddy  viscosity  form  of  the  standard  K  —  e  model  (see 
Launder  and  Spalding  1974)  with  a  coefficient  C*  that,  for  the  second-order  closures  consid¬ 
ered  herein,  is  close  to  the  traditional  value  of  0.09.  When  quadratic  terms  are  maintained, 
we  have 


(65) 


K,  =  -s-j  -  (smx  +  s-tw-ki )  +  2  (sjsfc  -  ls;n5,;n«ii) 

which,  in  an  inertial  frame,  yields 

th  =  Iks,,  -  2 c;—Sij  -  fa ~(s,kukJ  +  sjk&ki)  +  fa~  (s,kskj  -  \smnsm 

where  fa  =  \g2(2  —  C4)  —  C2)  and  fa  =  g2{2  —  C3)  —  C2).  This  is  identical  in  form 

to  the  nonlinear  K  —  e  model  of  Speziale  (1987)  when  convective  effects  are  neglected.  The 
nonlinear  K  —  e  model  is  given  by 


)  (66) 


K2— 


,  A'3  74 


1  4 


t„  =  -I<S,]-2C,~S,,-4CDC2—[S,J--SmmS,J 


-4 CDC2^  [sikSki  -  ^5 'mnSmnSi,) 


(67) 
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where  Co  and  C M  are  constants  and 


i  _  DS* 

ij~  Dt 


dvip 

dxk  k 3 


fain 

Q - 

ox  k 


is  the  Oldroyd  derivative  of  S.  When  DSij/Dt  is  neglected,  it  is  a  simple  matter  to  show 
that  (67)  is  of  the  same  form  as  (66). 

Eq.  (66)  is  also  of  the  same  form  as  the  anisotropic  eddy  viscosity  models  of  Yoshizawa 
(1984)  and  Rubinstein  and  Barton  (1990)  with  one  exception:  it  does  not  contain  a  term  of 
the  type 

f  Uik&kj  -  ■  (69) 

This  kind  of  term  was  eliminated  by  Speziale  (1987)  because  it  yields  an  erroneous  prediction 
for  isotropic  turbulence  subjected  to  a  solid  body  rotation.  More  precisely,  a  term  of  the 
form  (69)  predicts  that  an  initially  isotropic  turbulence  develops  Reynolds  stress  anisotropies 
when  subjected  to  a  solid  body  rotation  -  a  result  that  is  in  contradiction  of  physical  and 
numerical  experiments  (see  Wigeland  and  Nagib  1978  and  Speziale,  Mansour  and  Rogallo 
1987).  It  is  encouraging  that  this  consistency  constraint,  concerning  the  vanishing  of  (69) 
to  the  second  order  in  the  mean  velocity  gradients,  has  now  been  obtained  directly  by  a 
systematic  analysis  of  the  Reynolds  stress  transport  equation. 

It  is  interesting  to  note  that  the  quadratic  approximation  (66)  to  the  three-dimensional 
algebraic  stress  model  is  of  the  same  tensorial  form  as  the  full  nonlinear  two-dimensional 
model  (60).  However,  due  to  the  strain-dependent  coefficients  in  (60),  more  physics  is  ac¬ 
counted  for.  For  example,  Eq.  (60)  properly  predicts  that  in  a  rapidly  rotating  frame  -  where 
12,  and  hence  £,  tends  to  infinity  -  6*-  and  V  go  to  zero.  This  is  the  restabilization  effect 
of  turbulence  in  a  rapidly  rotating  frame  that  the  usual  anisotropic  eddy  viscosity  models 
are  unable  to  predict  (see  Speziale,  Gatski  and  Mac  Giolla  Mhuiris  1990).  Hence,  we  feel 
that  the  two-dimensional  algebraic  stress  model  (60)  should  be  tried  in  the  future  as  an 
alternative  to  the  more  commonly  used  quadratic  anisotropic  eddy  viscosity  models.  Some 
applications  of  (60)  will  be  considered  in  the  next  section. 


5.  ILLUSTRATIVE  EXAMPLES 

In  order  to  demonstrate  the  efficacy  of  the  regularized  algebraic  stress  model  (60)  derived 
in  this  study,  we  will  consider  a  few  applications  to  non-trivial,  two-dimensional  turbulent 
flows  involving  shear  and  rotation.  We  will  base  our  calculations  on  the  Speziale,  Sarkar  and 
Gatski  (1991)  second-order  closure  (hereafter,  referred  to  as  the  SSG  model).  The  regularized 


13 


algebraic  stress  model  can  be  rewritten  in  the  form 


6(1  +  t]2)aiK 


T,J  SKSii  3  +  if  +  6CV  +  6C2 


s-j  +  (S1WJ  +  s-„w^)  -  2  (s-ksij  -  isssj^)] 


(70) 

where  o;i  =  (c 2  —  /( C3  —  2)  and  S*,  W”,  7?  and  (  are  as  defined  in  (25),  (26),  and 
(57).  Some  remarks  are  needed  concerning  how  V/e  and  II are  evaluated  in  the  model 
coefficients  (21)  for  the  SSG  model.  V/e  is  calculated  using  (55)  which  is  formally  valid 
for  equilibrium  homogeneous  turbulent  flows;  we  take  //<,  «  0.11  which  is  the  universal 
equilibrium  value  predicted  by  the  SSG  model  for  two-dimensional  homogeneous  turbulence 
(see  Speziale,  Sarkar  and  Gatski  1991).  This  yields  the  following  choice  of  constants 


Ci  =  6.80,  C2  =  0.36,  C3  =  1.25,  C4  =  0.40,  g  =  0.233 


(71) 


for  the  explicit  algebraic  stress  model  corresponding  to  the  SSG  second-order  closure. 

The  explicit  algebraic  stress  model  (70)  is  solved  in  conjunction  with  modeled  transport 
equations  for  K  and  e.  These  equations  are  of  the  same  general  form  as  those  used  in  the 
K  —  e  model  and  are  given  by: 


where  vt  =  CMA'2/e,  =  0.09,  <jk  —  1,  <re  =  1.3,  Cc\  =  1.44  and  Ce2  =  1.83.  Hence,  the 
explicit  algebraic  stress  model  derived  herein  formally  constitutes  a  two-equation  turbulence 
model;  the  K  —  e  model  is  recovered  when  (70)  is  linearized  with  respect  to  the  mean  velocity 
gradients.  While  this  explicit  algebraic  stress  model  is  a  two-equation  model,  it  incorporates 
much  more  physics  than  the  K  —  e  model  since  it  is  consistent  with  second-order  closures  in 
the  limit  of  equilibrium  homogeneous  turbulent  flows. 

The  first  example  that  we  will  consider  is  the  case  of  homogeneous  shear  flow  in  a  rotating 
frame.  An  initially  decaying  isotropic  turbulence  is,  at  time  t  =  0,  subjected  to  a  uniform 
shear  rate  S  in  a  reference  frame  rotating  steadily  with  angular  velocity  0.  The  corresponding 
mean  velocity  gradient  tensor  is  given  by 


dvi 

dxj 


0  5  0  \ 
0  0  0 
0  0  0  / 


(74) 


and  the  angular  velocity  of  the  reference  frame  is  given  by  ft,  =  (0,0,  ft)  (see  Figure  1).  In 
Table  1,  the  equilibrium  values  predicted  by  the  new  explicit  algebraic  stress  model  (ASM) 
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for  homogeneous  shear  flow  are  compared  with  the  experimental  data  of  Tavoularis  and 
Corrsin  (1981)  as  well  as  with  the  predictions  of  the  SSG  model  and  the  standard  K  —  e 
model.  From  these  results  it  is  clear  that:  (a)  the  new  explicit  ASM  performs  far  better  than 
the  standard  K  —  e  model,  and  (b)  the  new  explicit  ASM  yields  results  in  close  proximity 
of  the  full  SSG  second-order  closure  model.  The  small  differences  between  the  new  ASM 
model  and  the  full  SSG  model  arises  from  two  sources:  the  neglect  of  the  quadratic  return 
to  isotropy  term  and  the  Pade  approximation  that  was  implemented  to  regularize  the  model. 
However,  it  is  clear  from  these  results  that  the  new  regularized  ASM  model  constitutes  an 
excellent  approximation  to  the  full  SSG  model  for  equilibrium  flows. 

Now,  we  will  compare  the  predictions  of  these  three  models  for  the  time-evolution  of  the 
turbulent  kinetic  energy  ( K *  =  K/Ko,  t*  =  St)  with  the  large-eddy  simulations  (LES)  of 
Bardina,  Ferziger  and  Reynolds  (1983)  for  rotating  homogeneous  shear  flow.  In  Figures  2(a)- 
2(c),  the  model  predictions  for  three  rotation  rates  (fI/5  =  0,  0/5  =  0.5  and  0/5  =  —0.5) 
are  compared  with  the  LES  which  is  for  an  initial  condition  of  e0/SK0  =  0.296.  From  these 
results  it  is  clear  that  the  new  algebraic  stress  model  does  an  excellent  job  in  capturing  the 
trends  of  the  LES.  It  yields  results  that  are  far  superior  to  the  standard  K  —  e  model  and 
are  in  close  proximity  to  the  SSG  model.  The  main  discernible  difference  between  the  new 
algebraic  stress  model  and  the  SSG  second-order  closure  is  during  the  early  transient  where 
it  responds  more  abruptly  to  the  application  of  the  shear  since  it  does  not  account  directly 
for  relaxation  effects.  However,  it  is  important  to  note  how  the  new  ASM  is  far  superior  to 
the  standard  K  —  e  model  in  responding  to  changes  in  the  rotation  rate.  This  is  illustrated  in 
Figures  3(a)-3(c)  where  it  can  be  seen  that  the  new  algebraic  stress  model  is  able  to  capture 
the  effect  of  rotations  on  homogeneous  shear  flow.  In  contrast  to  these  results,  the  standard 
K  —  e  erroneously  yields  the  same  results  for  all  rotation  rates  -  a  deficiency  tied  to  the 
Boussinesq  eddy  viscosity  hypothesis. 

The  last  example  that  we  will  consider  is  the  case  of  fully-developed  turbulent  channel 
flow  subjected  to  a  spanwise  rotation  with  constant  angular  velocity  0  (see  Figure  4).  We  will 
consider  the  experimental  test  case  of  Johnston,  Halleen  and  Lezius  (1972)  for  a  Reynolds 
number  Re  =  11,500  and  a  rotation  number  Ro  =  0.21  (here,  Re  =  UqH/v  and  Ro  = 
ilH/U0  where  Uq  is  the  bulk  mean  velocity).  In  Figure  5(a)  the  prediction  for  the  mean 
velocity  profile  obtained  from  the  SSG  second-order  closure  by  Speziale,  So  and  Younis 
(1992)  is  compared  with  the  experimental  data  of  Johnston,  Halleen  and  Lezius  (1972).  The 
computations  were  done  using  law  of  the  wall  boundary  conditions.  In  Figure  5(b),  the  same 
comparisons  are  made  for  the  new  explicit  algebraic  stress  model  (60)  derived  herein.  It  is 
clear  from  these  results  that  the  new  ASM  yields  an  asymmetric  mean  velocity  profile  that 
is  very  similar  to  that  obtained  from  the  full  second-order  closure.  Due  to  the  use  of  wall 
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functions  -  as  well  as  the  neglect  of  roll  instabilities  (see  Speziale,  So  and  Younis  1992)  -  the 
specific  quantitative  comparisons  are  not  that  good.  The  important  point,  however,  is  that 
the  new  algebraic  stress  model  yields  an  asymmetric  mean  velocity  profile  in  close  proximity 
to  that  obtained  from  a  full  second-order  closure.  On  the  other  hand,  simpler  two-equation 
models  like  the  standard  K  —  e  model  erroneously  predict  that  the  mean  velocity  profile  is 
unaffected  by  a  system  rotation  and  remains  symmetric  (see  Speziale  1991).  The  improved 
predictions  of  the  new  explicit  ASM  in  rotating  shear  flows  arises  from  the  incorporation  of 
rotational  strains  through  the  terms  W*  and  £. 

6.  CONCLUSION 

Explicit  algebraic  stress  models  for  two  and  three  dimensional  turbulent  flows,  in  non- 
inertial  frames,  have  been  obtained  for  a  hierarchy  of  second-order  closure  models  that  are 
tensorially  linear  in  the  Reynolds  stress  anisotropy.  These  models  were  obtained  using  the 
standard  local  equilibrium  hypothesis  and,  therefore,  constitute  the  explicit  solution  to  the 
traditional  algebraic  stress  models  generalized  to  include  non-inertial  effects  and  a  range  of 
pressure-strain  models.  They  also  formally  represent  the  equilibrium  states  predicted  by  this 
hierarchy  of  second-order  closures  in  homogeneous  turbulent  flows.  A  direct  examination  of 
these  explicit  models  has  shed  new  light  on  the  limitations  of  the  traditional  algebraic  stress 
models  in  applications  to  complex  turbulent  flows.  For  localized  strain  rates  that  are  large, 
traditional  algebraic  stress  models  can  become  singular.  These  models  need  to  be  regularized 
-  a  task  that  can  be  achieved  by  means  of  a  Pade  approximation  as  demonstrated  in  this 
paper. 

The  results  of  this  study  have  presented  the  first  definitive  evidence  as  to  why  the  tradi¬ 
tional  algebraic  stress  models  are  ill-behaved,  yielding  numerical  problems  in  many  applica¬ 
tions.  These  models  should  eventually  be  abandoned  in  favor  of  regularized  versions  of  the 
explicit  algebraic  stress  models  derived  herein.  The  full  three-dimensional  form  of  the  explicit 
models  is  rather  complicated  and  work  is  underway  to  simplify  them  by  means  of  a  ratio¬ 
nal  approximation  procedure.  However,  in  the  mean  time,  the  regularized  two-dimensional 
form  (60)  can  be  used  in  practical  applications.  This  simplified  model  is  formally  valid  for 
two-dimensional  equilibrium  flows  and  collapses  to  the  recently  proposed  anisotropic  eddy 
viscosity  models  in  the  limit  of  small  strain  rates.  However,  while  its  general  tensorial  form 
is  as  simple  as  the  commonly  used  anisotropic  eddy  viscosity  models,  it  incorporates  much 
more  physics  since  the  coefficients  depend  nonlinearly  on  both  rotational  and  irrotational 
strains.  A  version  of  (60)  based  on  the  SSG  second-order  closure  was  applied  to  rotating 
shear  flows  herein  with  encouraging  results.  These  developments,  when  combined  with  some 
recent  improvements  in  the  modeling  of  the  turbulent  dissipation  rate,  can  lead  to  a  new 
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generation  of  two-equation  models  that  can  serve  as  a  useful  companion  to  second-order 
closures  in  the  calculation  of  complex  turbulent  flows. 
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APPENDIX  A 


The  elements  of  the  matrices  H^\  and  J~,\  are  given  in  terms  of  t]i  as 
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The  elements  of  the  matrix  A  are  a  linear  combination  of  the  elements  of  the  matrices 
H-yx  and  JlX,  as  well  as  the  identity  matrix.  The  matrix  structure  is  then  given  by 
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with  the  nonzero  elements  determined  from 
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Equilibrium 

Values 

Standard 

K  —  e  Model 

New  ASM 
Model 

SSG  Model 

Experimental 

Data 

bu 

0 

0.204 

0.218 

0.20 

D 

-0.217 

-0.157 

-0.163 

-0.15 

^22 

0 

-0.149 

-0.146 

-0.14 

&33 

0 

-0.055 

-0.072 

-0.06 

SK/e 

4.82 

6.02 

5.76 

6.0 

Table  1.  Comparison  of  the  model  predictions  with  the  experimental  equilibrium  values  in 
homogeneous  shear  flow  measured  by  Tavoularis  and  Corrsin  (1981). 
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-  New  Explicit  ASM 

- SSG  Model 

- K—e  Model 


o  LES 


I'igure  2.  Time  evolution  of  the  turbulent  kinetic  energy  in  rotating  homogeneous  shear 
flow:  Comparison  of  the  model  predictions  with  the  large-eddy  simulations  of  Bardina  et  al. 
(1983).  (a)  n/S  =  0,  (b)  n/S  =  0.5  and  (c)  fl/S  =  -0.5. 
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—  New  Explicit  ASM 

-  SSG  Model 
--  K—e  Model 
o  LES 


Q/S  =  -0.5 


Figure  2.  Time  evolution  of  the  turbulent  kinetic  energy  in  rotating  homogeneous  shear 
flow:  Comparison  of  the  model  predictions  with  the  large-eddy  simulations  of  Bardina  ct  al. 
(1983).  (a)  f yS  =  0,  (b)  ft/5  =  0.5  and  (c)  ft/5  =  -0.5. 


(a)  LES 


Figure  3.  Time  evolution  of  the  turbulent  kinetic  energy  in  rotating  homogeneous  shear  flow: 
(a)  large-eddy  simulations  of  Bardina  et  al.  (1983),  (b)  new  explicit  algebraic  stress  model, 
and  (c)  standard  I(  —  e  model. 
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Figure  3.  Time  evolution  of  the  turbulent  kinetic  energy  in  rotating  homogeneous  shear  How: 
(a)  large-eddy  simulations  of  Bardina  el  al.  (1983),  (b)  new  explicit  algebraic  stress  model, 
and  (c)  standard  K  —  e  model. 


o 


Q 


Figure  4.  Schematic  of  fully-developed  turbulent  channel  flow  in  a  rotating  frame. 
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Figure  5.  Comparison  of  the  model  predictions  (-)  for  the  mean  velocity  profile  with  the  ex¬ 
perimental  data  (□)  of  Johnston  et  al.  (1972)  for  rotating  channel  flow  (Re  =  11,500,  Ro  = 
0.21):  (a)  SSG  second-order  closure  (taken  from  Speziale  et  al.  1992)  and  (b)  new  explicit 
algebraic  stress  model. 
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